Absence of Bose condensation in certain frustrated lattices 
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We study hard-core bosons on a class of frustrated lattices with the lowest Bloch band having a 
degenerate minimum along a closed contour in the reciprocal space. We suggest that the ground 
state of the system is given by non-condensed state, which may be viewed as a state of fermions 
subject to Chern-Simons gauge field. At fixed density of bosons, such a state exhibits domains of 
incompressible liquids. Their fixed densities are given by fractions of the reciprocal area enclosed 
by the minimal energy contour. 



There is a recent spike of interest in two dimensional 
lattices with flat bands both in the context of spin liq- 
uids [1-4] and cold atomic gases in optical lattices [H|-[9[ . 
The archetypical example is a Kagome lattice, where one 
of the bands (the lowest one for antiferromagnetic sign of 
hopping integral) is completely non-dispersive Ed. it 
was shown [12j that below certain critical filling fraction 
vq (e.g. = 1/9 for Kagome lattice) the groundstate of 
hard core bosons is infinitely compressible. The reason is 
in the exact degeneracy of non-interacting basis, allow- 
ing to form a certain number of non-overlapping localized 
states. Above the critical concentration the compressibil- 
ity is finite, but U(l) symmetry remains unbroken signi- 
fying the absence of Bose-Einstein condensation (BEC). 
In the language of s pin models, such a groundstate is 
known as spin-liquid 



13H22]. 



In this letter we show that BEC is absent in a very 
different family of lattice models. Namely we focus on 
2D lattices which lowest energy band exhibits degenerate 
minima along a closed contour in the reciprocal space, 
Fig.[TJ The simplest realizations of this scenario are given 
by lattices with two cites per unit cell. As an example 
consider a honeycomb lattice with nearest t\ and next- 
nearest ti hopping, Fig. El The common feature of all 
such lattices is the Hamiltonian of the form 



H = t±T + t 2 T 2 




(1) 



where the matrix structure is in A/B sublattice space 



and G 



* k e j with the three lattice 



^=1,2,3' 

vectors connecting a site of sublattice A with three 
nearest neighbors of sublattice B. Three unit vectors ej 
form 120 ° angles with each other, Fig. O The Hamil- 
tonians of the form (pQ) are not limited, though, to the 
honeycomb lattice. A generic oblique lattice with three 
distinct nearest and three distinct next-nearest hopping 
integrals is described by Eq. (pp), if two conditions are 
imposed on six hopping constants [23[ (variety of other 
lattices give rise to Hamiltonians of the form (pQ)). 

The two energy bands of the Hamiltonian (pQ) are given 
by £^ T) = TM|G k | + t 2 |G k | 2 . The lowest energy band 




FIG. 1. (Color online) Lowest energy band of the honeycomb 
lattice with t 2 /ti = (a), 0.3 (b), 0.5 (c), 0.9 (d). The mini- 
mal energy contour C is shown in light gray (blue). 



E^ exhibits a degenerate (flat) minimum along the con- 
tour C in the reciprocal space given by |GJ = |£i|/2£ 2 - 
For the honeycomb lattice this condition [24] is satisfied 
for t<i > \ti |/6, Fig. [TJ A similar dispersion relation ap- 
pears in the contex t of pa rticles with isotropic Rashba 



spin-orbit coupling [25 



The issue of Bose condensation for particles with such 
a dispersion relation is a non-trivial one. On the non- 
interacting level there is no transition at any finite tem- 
perature. This is due to the square root, [E — Ec) -1 ' 2 : , 
divergence of the single particle density of states (DOS) 
near the bottom of the band. Such behavior of DOS high- 
lights similarities with one-dimensional systems, where 
the groundstate of strongly repulsive bosons is given by 
the Tonks- Girardeau gas of free fermions [3ll435|. Here 
we show that the effective fermions picture describes the 
groundstate of hardcore bosons on our 2D lattices as well. 
The key observation [29| is that the chemical potential of 
free fermions with the dispersion relation of Fig. [T] scales 
as fiF cx v 2 at small enough filling factors v. This is 
energetically favorable in comparison with BEC in one 
of the states along the contour C. The chemical poten- 



2 




tial of the latter scales as [ib oc v. To build a fermionic 
state of Bose particles one uses Chern-Simons flux attach- 
ments familiar in the context of the fractional quantum 
Hall effect [36H38|. In the mean-field approximation this 
leads to fermions subject to a uniform magnetic flux 4ttu 
per unit cell superimposed with a particular staggered 
HaldanejZlj flux arrangement. Due to the Pauli prin- 
ciple, the fermions automatically incorporate hard-core 
condition and thus may be considered as non-interacting. 
This is the peculiarity of the lattice model with on-site 
interactions. In a continuum model with spin-orbit inter- 
actions, particles with opposite spins still interact 28l.l29|. 

Our main conclusions are as follows: the groundstate 
energy of hard-core bosons is minimized for a set of frac- 
tionally quantized filling fractions, determined by the area 
Ac (normalized to a total area of the first Brillouin zone) 
enclosed by the band minima contour C in the reciprocal 
plane |42| 



Ac 



21 + 1 + 7/tt ' 



Z = 0,1,... 



(2) 



where Berry phase is 7 = 0, if C encircles the T point of 
the Brillouin zone and 7 = 7r, if C encircles K and K' 
points, Fig. [TJ States with such filling factors are gapped 
and incompressible. For density in between the fraction- 
ally quantized fillings the system splits into domains with 
filling factors v\ and v\+\. Those states are U{1) sym- 
metric states of composite fermions (CF) occupying the 
lowest Landau level (or rather the lowest branch of Hof- 
stadter butterfly). The corresponding phase diagram is 
schematically depicted in Fig. 03 

Hamiltonian may be written in terms of bosonic cre- 
ation and annihilation operators frj, 6 r , with the hard- 
core condition = (b r ) 2 = 0, {6j,6 r } = 1 (where 
{., .} denotes anti-commutator), while operators on dif- 
ferent sites commute. For, e.g., honeycomb lattice the 



FIG. 3. (Color online) Phase diagram of hard-core bosons on 
a honeycomb lattice. CF and BEC are composite fermion and 
Bose condensate states respectively. Also shown incompress- 
ible states with fractionally quantized filling fractions vi. 



Hamiltonian takes the form: 



where the vectors ej and a^, j 



(3) 

1, 2, 3 are shown in 
Fig. O Chemical potential, /i, is related to the average 
on-site occupation v through an equation of state. 

Motivated by the observation that for divergent DOS 
the fermionic chemical potential is lower than that of 
the Bose condensate, we proceed with the Chern-Simons 
transformation 36l440| . To this end we write the bosonic 
operators as 



rVr arg[r-r jn r 



(4) 



where the summation runs over all sites of the lattice. 
Since the bosonic operators on different sites commute, 
the newly defined operators c r and c\ obey fermionic 
commutation relations. Also notice that the number op- 
erator is given by n r = c\c v . Upon transformation ((4]) 
hopping terms of the Hamiltonian (|3j) acquire phase fac- 
tors e 1 ^" ^r" y ,v n r" ? w here r " 5 r', r is a scanning angle of 
the link (r',r) seen from the lattice site r" . In terms of 
the fermionic operators the Hamiltonian (j3j) reads as 



V.J 



. r c r _|_ a . 



,r ,r+e; "'r / 



(5) 



JT I <Pr',r,r + a • n r' 



j ^ + H.c. 



Notice that the hard-core condition is taken care of 
by the Pauli principle and thus fermions may be con- 
sidered as non-interacting. A fermion hopping along 
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a closed loop £ on a lattice acquires a phase factor 



n 



<r'r">€£ 



= n r 



To 



analyze the consequences of these factors we adopt the 
mean-field approximation proven to be effective in the 
context of the fractional quantum Hall effect 36, l38| . 
n r ~ (n r ) = v. The phase factor takes the form 
e i<pc» = n r e^, where <p r = £ <r / r „ >GjC 0r,rT" is the 
total angle obtained by scanning the loop C from the 
point r. It is clear that cp r = if r is outside of the 
loop and cp r = 2tt if r is inside. If the reference point 
r is exactly on the loop/polygon, ip r is the angle of the 
polygon corresponding to vertex r. Therefore for a poly- 
gon C with s vertexes enclosing m lattice sites, the phase 
factor is cpc = (s — 2)tt + 2i\m. This formula has a sim- 
ple meaning: for any triangle (big or small) the flux is 

The unit cell of the honeycomb lattice consists of four 
triangles, Fig. [2j resulting in the average flux = Attv 
per unit cell. It is clear however that this flux is dis- 
tributed non uniformly: the half of the unit cell which 
contains site of sublattice B carries 37r^, while the other 
half 7tu. There is thus a modulation ±7rz^ between the 
two halves of the unit cell. It is convenient to divide 
this modulation between the three small triangles, each 
carrying additional flux (f>H = ttu/3 = <l>/12. Therefore 
the mean- field treatment of the Chern- Simons phase re- 
sults in non-interacting fermions subject to a constant 
magnetic field $ = Attv superimposed with the staggered 
Haldane phase 41] <pH = 3>/12. Below we analyze conse- 
quences of this mapping. 

We first notice that the flux attachment described 
above preserves the special form of the Hamiltonian ([1]) 
H = t\T + ^2^ 2 , where T is the operator acting in 
the space (c£,cf) and G = e ie j ( k+Ar ), where A r 
is the vector potential of the average magnetic field with 
<I> = Akv flux per unit cell. Notice that in the pres- 
ence of the vector potential operators G and & do not 
commute. Employing Baker- Campbell-Hausdorff[43j for- 
mula, one finds GG+ = J2j e ia ^ k + Ar ) + ^^ +H.c. where 
— 771 = 772 = ??3 = 1 and 4>h — 3>/12. This is exactly the 
next-nearest neighbor hopping over sublattice A in pres- 
ence of the constant magnetic field and Haldane modula- 
tion <j> H . Similarly G^G = £V e ^ r (k+A r )-^^ + R Q ig 

the hopping along sublattice B in the same setup. There- 
fore the Haldane modulation, naturally appearing from 
the mean-field treatment of the Chern-Simons field, is 
necessary to preserve the form (pQ) of the Hamiltonian (a 
constant magnetic field without the modulation does not 
admit representation JT}). 

This observation greatly simplifies finding the spec- 
trum by reducing the problem to diagonalization of the 
operator T. We first analyze it in the semiclassical ap- 
proximation, applicable if the minimal energy contour C 
encloses relatively small fraction of the Brillouin zone. 
To this end we notice that the spectrum of the opera- 
tor G = Gk+A in Eq. (pQ) can be found using Onsager's 




FIG. 4. Hofstadter energy spectrum vs. filling fraction v £ 
[0,1/2], for £2 = £i/4, i.e. C is around the F point. Notice 
that the bottom of the Hofstadter spectrum is flat, which is 
a consequence of the fact that all Landau levels, (|8|). exhibit 
minima at the same energy E — — £?/4£2. 



relation [44( for Bohr-Sommerfeld quantization of quasi- 
classical cyclotron orbits in a magnetic field. Denote by 
Gj($), I = 1,2,..., the eigenvalues of Gk+A- Semiclassi- 
cally Gi(<&) can be found by: (i) considering the constant 
energy contours |Gk| = G of the bare operator in the re- 
ciprocal k space, and (ii) identifying Gj($) with energy 
of contours C\ having a normalized reciprocal area A(C\) 
given by: 



A(d) 



2 2tt J 2ir ' 



(6) 



where is a magnetic flux through a unit cell of the 
lattice, and 7 is the Berry phase fjji l46| . Finally, the 
spectrum of the Hamiltonian (pQ), which describes the lat- 
tice subject to the constant magnetic flux <I> and Haldane 
modulation <pn = $/12, is found in terms of Gj($) as 



E l ($) = -t 1 Gi($)+t 2 [Gi($)Y 



(?) 



Landau levels fljp) are non-monotonic functions of flux, 
see inset in Fig. [5] They all (apart from I = if 7 = 7r) 
reach the minimum at G = £i/2£2, i.e. exactly at the 
band minimum, where the corresponding cyclotron orbit 
coincides with the minimal energy contour C. Recalling 
that <!> = 47r^, one obtains the set of the filling factors 
v\, Eq. ([2]), where the Landau levels reach the minimum. 
Since the Chern-Simons transformation ([!]) attaches ex- 
actly one flux quantum per particle, fermions c r fully 
fill the lowest Landau level (LLL) at any lattice filling 
v. Therefore in the mean-filed approximation the many- 
body groundstate energy follows LLL. 

As two examples we consider the cases where C is 
close to the Y point, i.e. £2 > ti/6, and C is close 
to K and K' \ i.e. £2 ^> ti/2. In the first case, ex- 
panding near k = 0, we find |Gk| 2 ~ 9(1 — k 2 /2) 



-3.8 




FIG. 5. (Color online) Bottom part of Fig. H Thick (red) line 
represents the ground state energy per particle obtained nu- 
merically from the Hofstadter energy spectrum, Eq. ([TO]). Ar- 
rows show fractionally quantized filling fractions ([2]). Dashed 
line is the microscopic chemical potential exhibiting jumps at 
the fractionally quantized filling fractions. Inset: Semiclassi- 
cal Landau levels © as functions of the filling fraction v. 



(where 1/k is measured in units of lattice spacing), re- 
sulting in the constant energy contours with the normal- 
ized area A(C t ) = (3>/3/4tt) (1 - Gf/9). Taking into 
account quantization Eq. (j6]) with 7 = 0, we obtain 

Gi 



±\ 9- 2 v / 3$(Z + l/2) and 



-M/9 



2v / 3$(/ + l/2)+t 2 (9-2>/3$(Z-r-l/2)). 

(8) 

This semiclassical spectrum is shown in the inset in 
Fig. [5] In the second case, ex paneling around K point 
we find Gk ~ 3|fc|/2 and 7 = 7r[4l(, leading to 



(9) 



To go beyond the semiclassical approximation we con- 
sider the Hofstadter problem on the lattice, including 
Haldane modulation. For a rational flux <1> = Airp/q (p 
and q are positive integers) diagonalization of the opera- 
tor T reduces to Harper equation, which can be analyzed 
numerically. For such fluxes the spectrum splits onto 
q non-overlapping subbands, labeled by m = 1,2, ...q. 
The corresponding spectrum E m ^(^Y Fig. HJ acquires 
the form of the Hofstadter butterfly [47j. Notice the flat- 
ness of the lower edge of the spectrum, which reflects 
the divergent DOS at this energy. Figure [5] amplifies the 
lowest part of the Hofstadter spectrum. Landau levels, 
closely following Eq. J5J), are visible at small filling frac- 
tions. 

The mean- field Chern- Simons treatment maps the 
Hamiltonian (|3]) onto the system of non-interacting 
fermions in the Hofstadter spectrum. Since fermions 



have filling factor v = p/q, the many-body groundstate 
is given by occupying p (out of q) lowest subbands. The 
ground state energy per particle is given by 



N/q 



(10) 



m—1 k 



where TV is number of lattice sites. In Fig. [5] we show 
groundstate energy calculated this way vs. filling frac- 
tion. For small filling fractions it closely follows the semi- 
classical lowest Landau levels (j5J), exhibiting the minima 
at the fractionally quantized filling fractions v\, Eq. (|2]). 
This leads to a macroscopic chemical potential, exhibit- 
ing staircase shape with the jumps at the fractionally 
quantized filling fractions ([2]), see Fig. [5] The flat regions 
of the staircase imply phase separation into domains with 
fillings v\ and v\+\. 

There is direct mapping between the considered hard- 
core boson system Eq. (j3]) and XY lattice spin model. 
Indeed, one may express bosonic creation and annihila- 
tion operators in terms of on-site spin 1/2 operators 
bi^ — »• <7jf and (2n r — 1) — > a*, where n r = b+b r is the 
on-site density operator. Then the Hamiltonian Eq. (j3j) 
acquires the form 



H 



r+e 7 - 



cr 1 a, 



r+a,- 



H.c. 



(11) 

The z-magnetic field plays the role of the chemical poten- 
tial /i = 2/i, which is related to the average magnetization 
m = v — 1/2 through the equation of state \i(y). In this 
way, BEC state with broken ?7(1) symmetry translates 
into magnetically ordered state, while U(l) symmetric 
CF state is interpreted as spin-liquid. Comparison of CF 
groundstate energy, Eq. (fTQ|) , with the energy of the mag- 
netically ordered state of the classical spin model [iil. l49| 
shows that at v = 1/2 (i.e. no magnetic field) the for- 
mer wins for 1/5 < ^Ai < 1/2. This suggests the spin 
liquid ground state, Fig. [3j of the corresponding XY spin 
model, corroborating with recent simulations 18, 50|, 5l| . 
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